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ABSTRACT 

Using Wigner transforms of Green functions, we discuss non-equilibrium generalizations 
of spectral functions and occupation numbers. We develop methods for computing time- 
dependent spectral functions in conformal field theories holographically dual to thin-shell 
AdS-Vaidya spacetimes. 
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1 Introduction 

The problem of thermalization of a quantum system that is prepared in a highly excited 
state is ubiquitous in physics. In contrast to situations where the equilibrium state of 
a quantum system is slightly perturbed and its return to equilibrium can be described 
by linear response theory, there is no general approach to the problem of thermalization 
starting far from equilibrium. In some cases, e. g. when the system is dilute and composed 
of weakly coupled quasiparticles, kinetic theory provides a general solution. However, many 
quantum systems of interest are either dense or strongly coupled, or both, and kinetic theory 
does not apply. 

A widely studied example is the strongly interacting matter produced in high energy 
nuclear collisions [ ]. When thermally equilibrated, this matter is known as a quark-gluon 
plasma, but the collision is thought to produce an initial state composed of densely packed, 
highly excited, but locally coherent gauge field domains, often referred to as a "glasma" [2]. 
The thermalization of this initial state is an open problem of intense phenomenological 
interest [ ]. 

The late time evolution of the energy-momentum tensor of a thermalizing quantum 
system is described by the effective low-energy theory of energy and momentum, fluid dy- 
namics. Chesler and Yaffe [ I] and Heller et al. [5,6], have shown that fluid dynamics becomes 
a good approximation rather quickly. But this "hydroization" of a quantum system does 
not imply full or even approximate thermalization of the system, because the anisotropy of 
the energy- momentum tensor may remain quite large, and there is no assurance that the 
expectation values of other observables are near their thermal values. 

Time evolution from near-equilibrium initial states in a strongly coupled gauge theory 
with a holographic AdS/CFT dual has been extensively studied (see the review [7]). More 
recently, focus has moved to thermalization from far off-equilibrium initial states. One in- 
teresting quantity to consider in such studies is the holographic entanglement entropy [8], 
which was generalized to time-dependent backgrounds in [9], in particular the collapse of 
a homogeneous massless shell in the AdS space-time, leading to black hole formation in 
the bulk, which corresponds to thermalization in the dual field theory [10]. Thermaliza- 
tion was studied in more detail in [11-15] following the time evolution of entanglement 
entropy, Wilson loops and equal-time correlation functions. The most nonlocal quantity, 
the entanglement entropy, reached thermal equilibrium last, but still very fast, on a time 
scale saturating the causality bound [16]. A complementary approach to thermalization 
was investigated in [17-19], who considered thermalization after an anisotropic deforma- 
tion of the boundary spacetime [20] , in particular how fluctuations created near the horizon 
and dissipation come to a balance to satisfy a generalized fluctuation-dissipation theorem. 
Isotropization and thermalization were also investigated in [21], by considering a large set 
of bulk initial conditions. Rotationally symmetric isotropization including radial flow was 
recently investigated in [22]. 

Here we extend this previous work by considering time-like correlation functions. This 
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allows us to generalize concepts usually studied in equilibrium to non-equilibrium states, 
and to compute them holographically. Two important quantities of this kind are the spectral 
density and the occupation number distribution of states. Time-dependent versions of these 
quantities are discussed and illustrated in a simple example in section 2. . 

Our holographic computations will be carried out in the thin-shell AdS3 Vaidya space- 
time, 

ds 2 = - [r 2 - 6{v)R 2 ] dv 2 + 2dvdr + r 2 dx 2 . (1.1) 

which models time evolution following the sudden injection of energy in the dual field 
theory. We will develop two approaches for computing two-point functions with timelike 
separations. The first approach, described in section 3, uses the geodesic approximation, 
which has been employed for computing two-point functions of spacelike separated heavy 
operators in [11,13,14,23]. There are significant new challenges in the present case - real 
timelike geodesies do not reach the boundary of AdS space (and therefore cannot be used 
to connect timelike separated boundary points). An alternative approach could have been 
to use analytic continuation from Euclidean signature, but AdS- Vaidya does not allow 
a standard continuation of this kind. Indeed, the absence of a conventional Euclidean 
continuation with which to define Green functions is a general challenge in the study of 
non-equilibrium systems. 

Here we propose novel approaches to overcoming these obstacles within the AdS/CFT 
correspondence: (i) by defining a new, non-standard Euclidean continuation, (ii) by using 
complex geodesies, (iii) by developing a prescription for splicing propagators across the 
infalling shell. In a saddle point limit appropriate to heavy operators, we obtain an analytic 
expression for equal-space correlators that agrees precisely between all three approaches. 
Method (iii) does not require the probe to be heavy, and hence in Sec. 4 we use it to obtain 
numerical results for retarded two-point functions in momentum space. 

A complementary approach to the holographic non-equilibrium spectral function and 
the fluctuation-dissipation relation has been developed in [24,25]. They consider perturba- 
tive expansions (derivative and amplitude) away from thermal equilibrium, working first at 
leading order [2 i] (where they also study non-equilibrium shift effects) and then extending 
to higher orders [25]. In particular [25] also derives a generalized (holographic) fluctuation- 
dissipation relation for Wigner transformed spectral function and statistical function (the 
anticommutator Green function) which holds for non-equilibrium states which are perturba- 
tively connected to thermal equilibrium. It would be interesting and important to compare 
the approach of [24, 25] to the methods and context of this paper. 

2 Time-dependent spectral functions and occupation numbers 

2.1 Spectral function and occupation numbers in equilibrium 

Consider a system prepared in a state described by a density matrix p. In the interaction 
picture or in the Heisenberg picture, the density matrix is time-independent while the 
operators evolve in time. In order to study how thermalization proceeds, we will look 
at the evolution of correlation functions. The relevant correlators are the time-ordered 
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(Feynman) correlation funct 

iG F (x,t) =tr( / 3T(O(x,t)O(0,0))), (2.1) 

and the retarded correlation function 

zG R (x, t) = 0(t)tr(p [0(x, t), O(0, 0)]) . (2.2) 

If the system is translation invariant in space and time, we can consider the Fourier trans- 
forms of the correlation functions Gp(k,u) and G-&(k.,u). In Fourier space, the retarded 
correlation function is complex valued, and its imaginary part defines the spectral function 1 

p(k, u) = -2 Im G R (k, to). (2.3) 

The imaginary part of the time-ordered correlator contains additional information, which 
depends on the mean occupation number density n(k, u) of the Fourier modes. When the 
system is at thermal equilibrium, the correlation functions are related by the Fluctuation- 
Dissipation Theorem (FDT) 

[1 + 2n(k, u)] p(k, u) = -2ImG F (k,w) , (2.4) 

with the occupation number equal to the Bose (or Fermi) distribution n(k, u) = [e u ' T= fl)~ 1 . 
The authors of [17] have emphasized the usefulness of (2.4) as a measure of the degree of 
thermalization in holographic theories, in the context of a model where fluctuations are 
created near the stretched horizon and transported to the boundary. In these works the 
fluctuation-dissipation theorem has also been expressed as a relation between the sym- 
metrized and antisymmetrized correlation functions 

G s (x, t) = tr(p {0(x, t), O(0, 0)}) , G A (x, t) = tr(p [0(x, t), O(0, 0)]) , (2.5) 

so that 

G s (k, u) = [1 + 2n(k, w)] G A (k, w) . (2.6) 
2.2 Time-dependent generalizations 

In order to determine the degree of thermalization we would like to define time-dependent 
generalizations of the occupation number density n(k, uj) and the spectral function p(k, u). 
By definition, a time-dependent system in not time translation invariant, so the Feyn- 
man and retarded correlators explicitly depend on two instants of time: G ? f i r(x 2 , t 2 ; Xi, ti). 
Space is still homogenous, so the correlators depend only on xi — x 2 which can be Fourier 
transformed to momentum k. But a simple transform from time to frequency space u is no 
longer meaningful as we seek a time- dependent notion of occupation number and spectral 
density. We therefore resort to a well-known alternative - the Wigner transform. 

1 Thermal spectral functions were first studied in the AdS/CFT correspondence for Af = 4 super- Yang- 
Mills theory in [26,27]. 
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First introduce an average time T and a relative time t by 



t ^ T ~\ h = T + l' ( 2 - 7 ) 

and define T-dependent Green functions by keeping the average time T fixed and Fourier 
transforming in the relative time t: 



/oo 
dte^Gn^t), (2.1 
-oo 



where we have suppressed the spatial coordinates. The T-dependent spectral function is 
then defined as 

p(T,u) = -21mG R (T,u). (2.9) 

In this paper, our goal is to study how the spectral function evolves as a function of the 
average time T. In order to focus on information near T, we may consider small intervals 
of relative time t around t = 0. To this end, we add a Gaussian window function into the 
integrand 2 : 

r°° 2 

G a (T,co) = dte^e-i? G(T,t) , (2.11) 



where a controls the width of the time window, and G is a correlator. In the presence of a 
finite time window the T-dependent spectral function becomes 

p a (T,u) = -2lmG a , R (T,u) . (2.12) 

We can now define a generalized, time-dependent notion of occupation number density 
by the relation 



n a (k,T,u) = - 



lmG a!F (k,T,u) _ ' 



(2.13) 



_Im G^ R (k, T,u) 

In an interacting field theory we expect the occupation number density of the non-equilibrium 
system to evolve towards thermality (e^ w — 1) _1 (for bosons), so that the system eventually 
satisfies the fluctuation-dissipation theorem and becomes fully thermalized. 



2 This is similar in spirit to the Husimi distribution, which is a Gaussian smearing of the Wigner distri- 
bution on phase space [28]. It is also related to a well known method in signal analysis. For a signal s{t) 
with time dependent frequency, it is convenient to add a Gaussian weight factor into the integrand with a 
small window about a mean time instant t', 

s a (t',w)= f dte iut e~ i ^ 1 s(t) . (2.10) 

Then by studying s a (t' , uj) as a function of t' one can track changes in frequency This procedure is known 
as the Gabor transform, and it is a special case of a more general class of wavelet transformations (with 
more customized window functions called wavelets). In our case, we are studying functions of two time 
variables, and we are primarily interested in tracking changes as the untransformed variable T changes, 
hence we are setting if = 0. It would be an additional possibility to perform a full Gabor transform, in 
which case varying t' would allow us to study information at different time separations, but that is beyond 
the scope of this paper. 
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We note that the notion of occupation number density is only meaningful for momenta 
and frequencies where the spectral function has support: if a mode does not exist in the 
spectrum, it is meaningless to ask whether it is occupied. In physical quantities, the oc- 
cupation number density will always be accompanied by the spectral function, so that no 
contributions originate from regions in momentum space where the spectral function van- 
ishes. This is reflected in (2.4). The issue does not arise in the conformal field theories 
considered later in the paper but the simple quantum mechanical example in the follow- 
ing section has a discrete spectrum and in that case expressions like (2.13) need to be 
interpreted with care. 

As explained in appendix A, we can now follow the analysis in [38] to define a time- 
dependent entropy density 



s(k, u, T) = (n(k, T, lo) + 1) In (n(k, T, w) + 1) - (n(k, T, u)) In (n(k, T, u)) . (2.15) 



This definition is a natural extension of the conventional microcanonical definition of von 
Neumann entropy of a system, and reduces to the latter in equilibrium. It should be con- 
trasted with more specialized notions of entropy such as entanglement and Renyi entropies, 
which have been the subject of many studies of holographic thermalization. 

2.3 A simple example: the quenched harmonic oscillator 

To illustrate the out-of-equilibrium notion of a time-dependent spectral function in (2.9), 
and in (2.12) when a finite time window is introduced, we first consider the harmonic 
oscillator with characteristic frequency ojq. 

As we will eventually be interested in applying our methods to quantum field theories, 
it is useful to think about the harmonic oscillator as a + 1 dimensional quantum field 
theory. This quantum field theory has a single one particle state (v\0) with energy uq. 
Multiparticle states then appear as the excited states (a^) n |0). We will study the single 
particle spectral function obtained from the two point function of x(t) and the occupation 
number of the single energy eigenstate. It should be also noted that all of our results in 
this section generalize trivially to free quantum field theories with the substitution 



This is because free field theories are just a set of decoupled harmonic oscillators labeled 
by a spatial momentum k, and because the observables we consider are related to the two 
point correlation functions, which are diagonal in fc-space. 

The retarded Green function for the oscillator in its ground state is 




(2.14) 



where 



uo Vk 2 + m 2 . 



(2.16) 



G R (t) 



sin(a;ot) 



6{t). 



(2.17) 
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Writing 



we find the Fourier transform 



0(f) 



i 

2^ 



-iuit 



du- 



u + iQ- 



G R (u) = / dt G R (t) = — 



1 



Using 



1 



P 



the spectral density is given by 

p(u) = -2lmG R (u) 



71 



U 



[6(u - (Jq) - S(u + UJq)} . 



(2.18) 



(2.19) 



(2.20) 



(2.21) 



Following the general discussion of Section 2.2, we now introduce a time window into 
the Green function, 

pOO 2 

iut „-^r n u\ ( 2 .22) 



G a , R (u) = / dte^e-^G R (t) . 

J — oo 

Using (2.17), (2.18) and the known Fourier transform of a Gaussian, we find 



G a , R (u) 



dJ 



1 



a 



'Itx 



exp 



a 



(u'—u+Uq) 2 I — exp 



(to'— UJ— UqY 



The spectral function is then 



and using (2.20), we find 



p a ((j) = -2lmG aiR (u)) 



(2.23) 
(2.24) 



P<r(w) 



7T 

7T 
(J 



a 



2,71 



exp 



a 



a 



'2tj 



a 



exp — — (u + u o y 



[5 a (uj - (Jo) -5 a (u + uj )] . 



(2.25) 



The Gaussian time window simply results in a Gaussian smearing of the delta function 
peaks in the spectral density. 

Next we study a quenched harmonic oscillator, with frequency oji for t < t and frequency 
(jf for t > to- The time of the quench can be taken to be to = without loss of generality. 
The retarded Green function G R {t 2l t\) now depends on both t\ and t 2 , not only on their 
difference. It is given by 



G R (t 2 , h 



-isin[wi(t 2 -ti)]0(ta-ti) if h,t 2 <0, 

-^siniufih-t^eih-h) if ti,t 2 >0, 

^: sin[o;jti] cos[a;/t 2 ] — ^ cosf^ti] sin[cu/t 2 ] if t 1 < < t 2 

if t 2 < < ti 



(2.26) 
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The matching conditions at to = are such that for t\ < both Gr and dGn/dt2 are 
continuous across t 2 = 0, and for t 2 > both Gr and dGn/dti are continuous across 
tx = 0. 

In terms of the average time and relative time coordinates, ti = T — t/2, t 2 = T + t/2, 
the retarded Green function reads 



G R (T,t) 



'/*)*(*) 



sin w 



-T)] cos^f +T)]0(t) 
-^sin[ W/ (|+T)]cos[^(|-r)]^(t) 



if 
if 

if 



T < 0, | < -T 
T > 0, | < T 



- < T < 1 

2 ^ ^ ^ 2 



(2.27) 



We can then study its Wigner transform (2.8) and obtain the time-dependent spectral 
function (2.9). The main features of p(T, u) can be uncovered by inspecting our expression 
(2.27) for G(T,t). Consider, for instance, the very early time regime, T — > — oo. In this 
limit, G(T, t) is given by the first line in (2.27), and p(T, u) reduces to the spectral function 
for a time- independent harmonic oscillator with frequency Ui. Similarly, in the late time 
limit T — > oo, we recover the spectral function of a harmonic oscillator with frequency Uf. 
The time-dependent spectral function for a quenched harmonic oscillator thus interpolates 
between those of harmonic oscillators with the initial and final frequencies. To get a first 
idea of the intermediate-time regime, we can take the average time to be the transition 
time, T = t = 0. In this case (2.27) reduces to 



G R (Q,t) 



1 



sm[( ' ; )t] + 



1 



sin[( 



u f 



9{t) . (2.28) 



and we see that at T = the retarded Green function executes beats between the initial 
and final frequencies in relative time. 

The Fourier transform of the retarded Green function with respect to relative time can 
be carried out explicitly for general values of the average time. Taking the imaginary part 
then gives the following time-dependent spectral function, 



p(T,u) 



sin [2T(ui-L)i)} 



sin [2T(u)+uij)] 



I J • J. ) {6(2u -u)t + u f ) - 5{2u + Ui + u f )) 

+7T 



if T<0 



i 



[S(2u — Ui — Uf) — 5{2u + Ui — Uf)) if T = 



(2.29) 



UJf 



sin [2T(u-ui f 



. [2T{uj+w f )] 



(uj-u> f )((2u]-uj f ) 2 -ujj) (uJ+UJ f )((2ul+Ul f ) 2 -LU 2 ) 



if T>0 



We observe that the time-dependent definition (2.9) allows for negative values of the spectral 
function for u > 0. This is to be compared with the equilibrium spectral function that is 
positive definite for positive u. 

Although it is not necessary, it may be useful to introduce time windows as in (2.11) 
in order to define generalized time-dependent spectral functions that only depend on the 
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behavior near T. This will have a smearing effect, for instance on the various delta functions 
appearing in our quenched harmonic oscillator example as can be seen in Fig. 1. Observe 
in particular that the spectral function can be made positive definite for u > for an 
appropriate value of the width a (Fig. 1 right). This is similar to the situation with 
conventional phase space Wigner density which attempts to define a notion of particle 
density in the phase space of a quantum theory. The Wigner density is sometimes negative, 
but the associate Husimi density, obtained by smearing with an appropriate Gaussian is 
everywhere positive [28]. 




0) 



Figure 1: Spectral function of a quenched harmonic oscillator (ui = 2, Uf = 1) with 
Gaussian time window p a (T,u) at T = for a = 2 (left) and a ~ 1.48 (right). 

Examining the smeared spectral function for different times (Fig. 2) shows that p a 
evolves between the spectral functions of harmonic oscillators with the correct initial and 
final frequencies. 

Finally, we can compute a time-dependent occupation number, as defined in (2.13). The 
result for specific parameter values is shown in Figure 3. The occupation number quickly 
evolves to its asymptotic value, which is given by the expectation value of the harmonic 
oscillator number operator (for oscillator frequency Uf) evaluated in the ground state of an 
oscillator of frequency u^, which is the state of the system at the time of the quench. The 
oscillations around the final value arise from the smearing due to the finite time window - 
their size depending on the choice of the parameter a. 

3 The geodesic approximation in the AdS/CFT correspondence 

To use the methods developed above to discuss time dependent spectral functions in the 
AdS / CFT correspondence we need to first compute two point functions in the latter setting. 
In this section, we use a geodesic approximation to compute two-point functions with 
timelike separation in a two dimensional CFT dual to the AdS3-Vaidya spacetime (1.1). 
This approximation is accurate for heavy operators in the CFT. We will find simple analytic 
results for the special case of equal-space two-point functions, but our methods generalize 
to generic separations. These results allow us to compute an integral over spatial momenta 
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Figure 2: Time evolution of p a (T,u) (with uji — 2, Uf — 1) for o = 7 and T = 
— 12, —7.2, —2.4, 2.4, 7.2, 12 from left to right. At early (late) times the spectral function is 
the one of a time-independent harmonic oscillator with frequency Ui (u)f). 



of the time-dependent spectral function of high-dimension operators, as well as a weighted 
integral of time-dependent occupation numbers. 

Consider a scalar operator 0(x,t) with conformal dimension A in the dual CFT. The 
time-ordered two-point function (0(x±, ti)0(x2, t 2 )) is given by a path integral over all 
paths V that connect the two insertion points (a?i,ti) and (£2, £2) on ^ ne boundary: 

VVe- AL ^ with L(P) = J^ giw — — dX. (3.1) 

For large A, we can use a saddle point approximation, in which the sum over all paths can 
be approximated by a sum over all geodesies connecting the boundary endpoints [ 

(0(x 1 ,t l )0(x 2 ,t 2 )) ~ e ~ A£ ' ( 3 - 2 ) 

geodesies 

where £ denotes the geodesic length. However, due to contributions near the asymptotically 
AdS boundary, the geodesic length between two boundary points contains a universal diver- 
gence and needs to be renormalized. Throughout this Section, we will define a renormalized 
length 5C = C — 2 ln(ro), in terms of the bulk cut-off ro, by removing the divergent part of 
the geodesic length in pure AdS (see Eq. (B.16)-(B.17) in Appendix B). The renormalized 
two point function can then be approximated by 

(0(x 1 ,t 1 )0(x 2 ,t 2 )) Icn ^ e - ASC , (3.3) 

where 8L is the renormalised length of the geodesic that connects the points (xi,ti) and 
(x 2 ,t 2 ) on the boundary. 



10 



0.20 r 



0.15 - 




Figure 3: The plot shows the occupation number for frequency w = wj = 1 as a function 
of average time (here T is multiplied by a factor of 2 compared with the convention in the 
text) for a quenched harmonic oscillator, with final frequency Uf — 1 and initial frequency 
uii = 2, with smearing parameter a = \/8. The limiting values agrees with a conventional 
computation of the final occupation number in terms of a Bogoliubov transformation of the 
initial vacuum. 



The geodesic approximation works most straightforwardly in Euclidean spacetimes, or in 
static spacetimes that can be Wick rotated to Euclidean signature. In [13, 14], this approach 
was used in a time- dependent setting, to determine equal-time two point functions in the 
field theory dual of the Vaidya spacetime. This work was further extended in [ ] to non- 
equal-time two-point functions with spacelike separations. In all these cases, the relevant 
bulk geodesies V were spacelike, and the geodesic length C(V) was real. 

New complications occur when the points on the boundary are timelike separated, in 
which case we might expect to connect them using timelike geodesies. The main problem 
is that (real) timelike geodesies in an asymptotically AdS spacetime never extend to the 
boundary. For example, in the case of empty AdS3 with metric 

ds 2 = -r 2 dt 2 + ^- + r 2 dx 2 , (3.4) 
the timelike geodesies (e 2 > j 2 ): 

r(A) = \/ e 2 — j 2 cos(A) 

x(A) = x + tan(A) (3.5) 
t(X) = t + tan(A) 

satisfy r(A) < \/ e 2 — j 2 , VA. Thus they do not reach the boundary, where r — > oo. 

In fact, this is a specific realization of the general difficulty of defining correlation func- 
tions in time-dependent spacetimes that do not have a Euclidean continuation. Given that 
a general prescription is not available, we will here propose and investigate three potential 
approaches: 
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1. For time-independent spacetimes the standard prescription for vacuum correlators 
is to calculate in Euclidean signature and then analytically continue to Lorentzian 
signature. Since there is no obvious Wick rotation of the AdS-Vaidya spacetime, we 
will propose and use a non-standard Euclidean continuation in Sec. 3.1. 

2. We will consider complexified geodesies in Sec. 3.2. In this approach, geodesies be- 
tween timelike separated boundary points make excursions in a complexified AdS- 
Vaidya geometry before returning to real points on the boundary. 

3. In Sec. 3.3 we will directly connect the AdS and black hole retarded bulk-to-boundary 
propagators across the infalling shell, and will apply a saddle point approximation to 
the resulting integral expression. This approximation does not directly rely on the 
geodesic approximation but shares the same regime of validity of large scalar masses. 

For equal-space correlators, we will find the same simple analytic expression in all three 
approaches. 

3.1 Continuing to Euclidean signature 

In appendix B, we review how AdS and BTZ can be Wick rotated to Euclidean signature, 
and how various two-point correlation functions of high-dimension operators in the dual 
field theory can be computed using bulk geodesies. Here we develop a similar approach for 
AdS-Vaidya. 

3.1.1 Wick rotation of AdS-Vaidya 

Recall that the thin-shell AdS-Vaidya metric 

ds 2 = - [r 2 - 6{y)R 2 ] dv 2 + 2dvdr + r 2 dx 2 (3.6) 

represents an infalling shock wave of 'null' dust in AdS 3 . Outside the shock wave (v > 0), 
the metric is given by the BTZ solution 

ds 2 = -(r 2 - R 2 )dt 2 + + r 2 dx 2 with t = v - i_ m , ( 3 .7) 

while inside (v < 0), we recover the pure AdS3 metric 

dr 2 1 
ds 2 = -r 2 dt 2 + — + r 2 dx 2 with t = v + -. (3.8) 

While for static spacetimes, such as AdS3 or BTZ, a Wick rotation of the metric is straight- 
forward, this is no longer the case for an explicitly time dependent geometry. Thus we 
propose the following procedure: 

• We formally take the Vaidya metric to be a limit of a family of "spacelike Vaidya" 
metrics. The latter need not be solutions to the equations of motion with conventional 
matter - they are auxiliary objects used to regularize the computation. 
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• This spacelike Vaidya metric can be Wick rotated to Euclidean signature by a simul- 
taneous rotation of the time coordinate, as well as of an additional parameter. 

• In this auxiliary Euclidean spacetime, correlation functions can be computed in a 
geodesic approximation. The result can then be Wick rotated back. Taking the limit 
we find our proposed result for the conventional AdS- Vaidya space-time. 

While the results we find will look reasonable, it is not a priori obvious that this procedure 
should be valid. This is why, after discussing the present method in detail, we will also 
develop two alternative prescriptions. All three prescriptions match, suggesting that the 
method of considering a time-dependent space-time as a limit of others with standard 
Euclidean continuations may be more generally useful. 



Lorentzian 'spacelike' Vaidya 



Let E > R > and consider the 'spaceike' Vaidya metric 

2E dvdr dr 2 



ds 2 = - \r 2 - 6(v/E)R 2 } dv 2 + , + -= — /rnNn0 + r 2 dx 2 . 

1 K ' ' J ^/r 2 + E 2 -6(v/E)R 2 r 2 + E 2 - 6{v/E)R 2 

(3.9) 

Through the coordinate transformation 

t - ^arctanh ( §Jl + ^-=^1 J + ^arctanh (f ) for v > 

t + ^ll-y/l + fj foiV<0 

we recover the BTZ metric (3.7) for v > and AdS3 (3.8) for v < 0. The coordinate 
transformations are such that v — describes the geodesic of an infalling shell of 'spacelike' 
particles in the BTZ and AdS3 spacetimes. Observe that for E > 0, 9{y/E) = 9(v), so that 
in the E — > oo limit (3.9) reduces to the Lorentzian 'null' Vaidya metric (3.6). 



Euclidean Vaidya 



On the Lorentzian 'spacelike' Vaidya metric (3.9), we can now perform an analytic 
continuation, on the time coordinate z = iv, as well as on the parameter Z = IE. Without 
loss of generality, we can take Z > and find the metric 

d s 2 = W 2 - 9(z/Z)R 2 } dz 2 2ZdzdT dr 2 2&2 ( n) 

i J /„ 2 _ ^2 _ e(z/Z)R 2 r 2 - Z 2 -6(z/Z)R 2 v ; 



where as before 9(z/Z) = 0(z). Note that the radial coordinate r now runs from \J Z 2 + R 2 
to oo. 
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Letting 



the metric becomes 



f(rf = { f '° ( f ]2 = r \~ Z ' *<° (3.12) 

fout{r) =r 2 -Z 2 -R l z > 



ds 2 = firfdz 2 + ( Zdz - +r 2 dx 2 , (3.13) 

V f( r )J 

which is manifestly positive. 

3.1.2 Geodesic length in Euclidean Vaidya 

We now compute the length of geodesies that start at a point (r ,Xi,?/i) and end at 
(r , X2, 1/2), where r denotes the location of the regularized AdS boundary. We assume 
Ay = i)2 — yi > and, for computational simplicity, mostly focus on geodesies with 
Ax = X2 — xi = 0. Geodesies in the Euclidean Vaidya geometry can be constructed 
by gluing together at the shell, with the appropriate refraction condition, geodesies in Eu- 
clidean AdS and Euclidean BTZ. The latter have been worked out explicitly in Appendix B 
and read 

1 /r 

x±(r) = xo± — arctanh 

r \ r 

1 

z±{r) = yo ± — arctan 

1 (R 

-arctan ^- 

A-t(r) = A ± arccosh ( \\ ^ ^2" I (3-14) 

for z > 0, and 




_ . . 1 . , / S 2 

x± (r) = x ± — cos((x)y 1 - — , 



1 / S 2 I 1 / z 2 

z±{r) = y ±-sm(a) 



E v ' V r 2 Z Z V r 2 



A-t(r) = A ± arccosh ^—j (3.15) 

for z < 0. Here ± denote two separate branches (which appear above and below the turning 
points of AdS or BTZ geodesies) and xo, T±, y , X , x , S, cr, y , X are constants. 
We can distinguish three cases. 
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• If both endpoints occur before shell injection (y\ < y 2 < 0), there is a geodesic 
connecting them which is entirely within Euclidean AdS. Any additional geodesies 
would require the existence of geodesies in BTZ that connect two points on the shell, 
and these do not exist. 

• If both endpoints occur after shell injection (0 < y\ < j/ 2 )) 5 there is a geodesic 
connecting them which is entirely within Euclidean BTZ. 3 

• Finally consider geodesies with y\ < and y% > 0, and such that Ax = 0. They cross 
the shock wave once and hence we need to determine how they refract at the shell 
location in order to compute their length. 



Boundary conditions at the shell 



Consider a point P« n (P ou t) just inside (outside) the shell, and let the coordinate differ- 
ence between P; n and P ou t be AX = (Az, Ar, At). Take another point M on the shell z = 
and let the coordinate difference between P in (P out ) and M be AX in = (Az in , Ar in , At. 
(AX out = (Az out , Ar out , At out )), so that AX in + AX out = AX. Then the distance from P { 
to Pout via M is 



in I 



As = JAsI + ^/A? 



out 



7^A4 + (zAz in - ^) + r^Axl 

+]j flAzLt + (zA Zout - ^j^^j 2 + r%Ax - Ax in f , (3.18) 
where fi n and f ou t where defined in (3.12). We want to find the point M that minimizes 



3 We can check that no new geodesies are admitted in the presence of the shell as follows. Consider a 
geodesic with both endpoints in Euclidean BTZ, i.e. after shell injection (0 < y\ < 2/2)- If the geodesic has 
to extend into the z < region, it will have to cross the shell twice. It will therefore consist of two z + (r) 
branches, which are connected by a geodesic in AdS3. z~(r) thus needs to have a minimum, located at 

ZScos(cr) 

r = - (3.16) 



Z 2 - £ 2 sin 2 (a) 

For the special case of geodesies that connect points on the boundary with Ax — 0, one has cos(cr) = 0, 
such that tq = 0. However, since r > \J Z 2 + R 2 , z_ (r) reaches no minimum and the geodesic is thus 
completely in Euclidean BTZ. The renormalized geodesic length is worked out in Eq. (B.18) and reads 



SCbbtz — hi 



(3.17) 
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As. Extremizing this with respect to Arj n and Axi n we find 



[ ZA 



"'it. 



Ar r 

fin / fin 



As? 



1 f 7 , Ar on A 1 

/ 2 I ^^out 7 I 1 

y As oui \ J out J J out 

V^Lt 



(3.19) 



and therefore 



j-2 I ^ z in 

AXjr 



1 



1 An 



fin V AXj n AX{ n 



out 



1 Ar n , 



out 



. Az on4 

A.T ou ^ /out Ax ou t 



+ z 



. Azi, 



1 An 



2 / A^ OMf \ 2 + f g Azput 1 Ar out \ 

\Ax out J \ Ax out f out Ax out J 



In the AX — > limit, we obtain the refraction law 



£2 dx 
J dr 



Zf 



dz 
dr 



£2 dx 
J dr 



dx\ 2 



dz 
dr 



+ [Zf 



dz 
dr 



Zf 



dz 
dr 



fout Ax OU ( J 

(3.20) 



(3.21) 



out 



dx\ 2 



{ft) 



4 ( dz\* ( „ „dz 



^ \dr J ~^ \ ^ dr 



out 



that is supplemented by the continuity conditions x in (r*) = x out (r*) and z in (r#) = z out (r#) 
0, where r* is the value of the radial coordinate at which the geodesic reaches the shell. 



Geodesies crossing the shell 

If yi < < i/2, the geodesic starts in Euclidean AdS, crosses the shell once and ends in 
Euclidean BTZ. The geodesies are given respectively by the expressions (3.15) and (3.14). 
Imposing the boundary conditions results in the following relations 



%2 
V2 

yi 



x + -^arctanh ^7 



Do + -narctan 



x 

m - 



(R 2 -r 2 _) 



COS^cr J 
sin(a) 



Together with the refraction condition (3.21) that leads to 
i? 2 £ 2 cosV) = T 2 + T 2 _ 



ZTj sin(cr) 



Z 2 



— \ T 



R 2 ) 



Zr^(Tl-R 2 )(R 2 -Tl] 
R^Jrl-{Z 2 + R 2 ) 



(3.22) 



(3.23) 



n)(r 2 -r 2 
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and the continuity conditions at the shell 



1 , s I £ 2 1 , /T_ r 2 -T 2 

x + - cos(a)W 1 - — = x + -arctanh —J - _ , , 



1 . , / £ 2 1 1 / Z 2 
Vo+ — sm(er) W 1 — — + — W 1 - = 



yo + 



E 
1 



r 2 Z Z 



(3.24) 



(3.25) 



arctau 1 r 2 - i?2 V r 2_ F 2 ) + arctan y z 



* 'l-^±^i|-axctan 



i? 

z 

(3.26) 

these fix all unknowns (xo, J/o, ^o, l/o, cr, T + , T_ and r*) and determine the renormalized 
geodesic length 



= 



5C 



ESheii = hm [A+(r ) - A_(r ) -21n(r )l 



4 



, rl - T% \ I i £2 \ ! / 4 \ i 
—arctanh | \ / — — | + arctanh | ^ / 1 | + — In | — | + — In 



z 2 J 2 V r +-r 2 

(3.27) 

where in the last line we have used the continuity condition A + (r*) = A+(r*). 

When Ax = 0: x = x = X\ = X2, cr = n/2 and T_ = 0, the geodesic length (3.27) 



becomes 



1 I r 

^^EShcll = - In 



cos( 7 (r,)) y/r 2 -R 2 4 [(r,M(r,)) 2 + if 



2 \n + cos( 7 (r <! ))V^ 2 -^ 2r * [(r 2 - J R 2 )sin 2 ( 7 (r,)) + J R2] 
where r* is defined by the parametric relation 



(3.28) 



1 



(r*M(r,)) 2 + 1 
and 



2Zr,M(r,; 



r 2 _ £2 



- (r*M(r*)) 2 + 1 



Zsin(7(r ;> )) 



with 



M(r*) 
sin(7(r*)) 

1 



1 1 



Z 2 



-yi + z-z' 11 -- 

Rsm(Rtt(r*)) 
r* — cos(M](r*)) a/V 2 — i? 2 



^0"*) = 2/2 + ^arctan 



i* / (Z 2 + J R 2 )\ 1 

1 — —arctan 



R 



'R 
~Z 



cos(7(r*)) 
(3.29) 

(3.30) 
(3.31) 

(3.32) 
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Performing the double Wick rotation y\ = it\, y 2 = it 2 and Z = iE of (3.28) and taking 
the E — > oo limit, in order to recover the 'lightlike' Vaidya result, we finally find 



5C 



Shell 



In 



2 . (Rt 2 

— sinn 

R V 2 



("•OSh ( — 1 ti 



(3.33) 



3.1.3 Two-point functions of thermalising CFT 

As reviewed in appendix B for the cases of AdS and BTZ, by Wick rotating the Euclidean 
two point function Ge(x 2 , Vi', Xi, Vi), one finds the time-ordered (Feynman) two-point func- 
tion 4 

iG F (x 2 , t 2 ; Xi,ti) = G E (x 2 , it 2 ; x 1: it x ) . (3.34) 
The retarded two-point function can be obtained from the time-ordered one by 

G R (x 2 , t% xx, ti) = 9(t 2 - ti) [G F (x 2 ,t 2 ; xx, ti) + (G F (x 2 , t 2 ; £i, ti))*] . (3.35) 

From (3.33) we then find (the geodesic approximation of) the Feynman and retarded 
two-point functions in the thermalizing CFT dual to 3-dimensional AdS- Vaidya: 

g— iirA \ / g— iirA 



u-u 



iG,(U..r.t ,..,■) [ - w }e(-t 2 )9(-t 1 )+ | - - - — 7TV ^ \(-)(t->)(-)(t i 

e" 



sinh(f(t 2 -ti))| 



-i7rA 

+ I ~ 9A I d(h)9(-ti 

1 '| sinh (i|2)_ cosh (i|l) tl | 2A 



e -inA 



"' -sinh^-cosh^t^ 1 ^ ' (3 - 36) 
G R (t 2 ,x; tl ,x) = -2sH«mt2-h) ( V 2 ^ + md{t 



i 2 -ti| 2A ||sinh(f(t 2 -ti))| 2A 
0{t 2 )6{-h) 
sinh(^) -cosh(^)ti| 2 



(3.37) 



R 



In the limit R — > or t 2 — > + , Gn(t 2 ,x;ti,x) reduces to the retarded vacuum two-point 
function (B.27) with At = t 2 — t\ or At = —t\ respectively, while for t% — > 0™ to the thermal 
correlator (B.28) with At = t 2 . The above formulas are valid for non-integer values of A. 
For a discussion of integer values of A, see Appendix C. 

To compute time-dependent spectral functions and occupation numbers, we would need 
to Fourier transform the relative times and positions in the retarded and time-ordered 
two-point functions. This requires knowledge of the position space two-point functions for 
arbitrary separations, which are more cumbersome to obtain than the equal-space two-point 
functions we have explicitly computed so far. We will therefore focus on the information 



4 The e prescriptions are given in Appendix C.2. 



18 



we can extract from the equal-space two-point functions, i.e. the time-dependent spectral 
function integrated over all spatial momenta. 

Following the discussion in Section 2, we introduce relative and average time coordinates 
t = t 2 — ti and T = (ti + 1 2 )/2 and Fourier transform Gn(t2,x; ti, x) with respect to t. The 
imaginary part of the result is 

/+oo / p+oo \ 

dk p(k, T, u) = -4tt Im ( / dt e iut G R (x = 0, T, t) j , (3.38) 

where the additional factor of 2-7T with respect to the definition (2.9) appears because the 
retarded two point function is here given in position space, rather than momentum space. 

Numerical results for the spectral function integrated over momenta are plotted in Fig. 4 
for different times. As time evolves, the curves interpolate between the vacuum result 
(in dotted blue in panel (A)) and the thermal one (in continuous red). As the spectral 
function is odd in u, we only plot the range u > 0. Panel (B) plots, as a function of time 
and for different values of scaling dimension A, the slope (3 of the spectral function in the 
range of linear growth close to u = 0. For a fixed frequency u>, p(co,T) does not increase 
monotonically from the vacuum to the thermal value, but it oscillates. 




A=1.75 
A=2.25 
A=2.75 



(A) (B) 

Figure 4: Panel (A) plots the numerical results for the spectral function integrated over 
momenta p(T,uj) for different values of the mean time T (with A = 9/4, R — 15). When 
T — 1/R, the curve reduces to the the vacuum result (in blue dotted), while for T 3> 1/R, 
to the thermalised one (in red continuous). (The reference case T = is also shown in 
purple.) Close to the origin, the spectral function grows linearly in frequency. The slopes 
(3 of this regime of linear growth are plotted in panel (B) as a function of time and for 
various values of A. Similar results are obtained for larger values of A. The oscillations 
in the slope characterize how the spectral function interpolates between the vacuum and 
thermal results. 



As in the example of the quenched harmonic oscillator we discussed in Sect. 2.3, the time- 
dependent spectral function we have defined can be negative for positive u. Introducing 
an appropriate time window as in Eq. (2.11)-(2.12), we can smear the curves and make the 
spectral function positive definite for u > (see Fig. 5). 

From the time-ordered two-point function, using definition (2.13), we can also obtain 
time-dependent occupation numbers weighted by the spectral function and integrated over 
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(A) 



Ar(T) 

1.5 . 

\r 

1.0 ' 
0.| 



-1.0 -0.75 -0.5 -0.25 



(B) 



A=1.75 
A=2.25 
A=2.75 



Figure 5: Spectral function integrated over momenta (A) and slopes of the regime of linear 
growth (B) as in Fig. 4, with time window set by a = 2.9/ R. 



all spatial momenta n(T,u). These are plotted in Fig. 6 as a function of time, without (in 
(A)) and with (in (B)) the Gaussian smearing. 




(A) 




(B) 



Figure 6: (A) The integrated occupation number weighted with the spectral function. The 
times for which this quantity becomes negative corresponds to the same times for which 
the spectral function has a region that satisfies up(u, T) < 0. (B) The same quantity with 
a time-window set by o = 2.9 /R. 



3.2 Complexified geodesies 

To investigate the validity of the non standard Euclidean continuation used in Sect. 3.1, 
we here repeat the computation of the geodesic length in the 3-dimensional AdS-Vaidya 
spacetime using a different approach. While real timelike geodesies can not connect timelike 
separated points on the boundary, complex geodesies can. The latter have been studied 
in [30-32] and we review their construction in both AdS3 and BTZ in Appendix D. As we 
show below, the geodesic length obtained with this second method coincides with (3.33). 

As before, we can construct Vaidya geodesies that cross the shockwave by joining to- 
gether AdS and BTZ geodesies, with an appropriate refraction condition at the shell. 
For equal space geodesies, the relevant formulae are given in Eq. (D.5) and (D.16) with 
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j — J — 0. For v > 0, the coordinate transformation v = t — -^arccoth (jy leads to 



r 2 



(A) = VE 2 -R 2 sinh(A - A 2 + i/3 2 ) 



v 2 (X) 



t 2 — — arccoth 
R 



t 2 — —arccoth 
R 



E 
R 



tanh(A — A 2 + i(3 2 ) 



E + VE 2 -R\ /l. x . 
— tanh I -(A - A 2 + %p 2 , 



while for t> < 0, v is defined as v = t — \ and the geodesies are given by 

ri (A) = e sinh(A — Ai + 

£i(A) = n cothfA - Ai + i/3i) 

e 



«x(A) 



7~i coth 

e 



-(A - Ai + 



The refraction condition at the shell has been computed in [13, 14] and reads 



r'(A*) 


r'(A*) 


i? 2 


u'(A*) 


in V '(K) 


~~ 9 

out ^ 


x'(A*) 


x'(A*) 


= 0, 

out 


i/(A„) 


in 



(3.39) 
(3.40) 

(3.41) 



(3.42) 
(3.43) 

(3.44) 



(3.45) 
(3.46) 



where the second condition is trivially satisfied for equal-space geodesies. Continuity at the 
shell fixes 



esinh(A!) =VE 2 - R 2 sinh(A 2 



r 2 =— arccoth 
R 



E + VE 2 - R 2 
R 



tanh 



2 



T\ =- coth 

e 



Ai 
2 



(3.47) 
(3.48) 

(3.49) 



where we have defined Ai = A* — Ai + i$\ and A 2 = A* — A 2 + i(3 2 . The boundary conditions 



arc 



t 2 = t 2 arccoth 

R 

1 

ti = n + - 



E + VE 2 - R 2 
R 



(3.50) 
(3.51) 



and using the continuity condition (3.47), the refraction condition can be written as 
^ = EVE 2 - R 2 cosh(A 2 ) - e 2 cosh(A 1 ) + e 2 - E 2 + R 2 . 



(3.52) 
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Finally, we find the solution 

R [cosh (4^) - fisinh(^)] 



E 



m 

2 



and 



# cosh (4^) 



m 

2 



icosh (^) - sinh (^) 

and exp(— A 2 ) = sinh 



exp(Ai) = — — coth 



i cosh(^) -sinh (4^) 
Rt 2 



(3.53) 




-1, (3.54) 



such that 

<5£shcii 



1, f4\ 1 
In I — ) + - In 



2 
In 



4 \ 1 / 16e 2Aa e- 2Al 

2 ~ V^ 2 - R 2 J + ~ ~ 2 11 Ve 2 (^ 2 -^ 2 ) 

| sinh (^)-cosh^ t x 



(3.55) 



which corresponds exactly to the result (3.33) that we found using the analytic continuation 
to Euclidean signature. 

3.3 A saddle point approximation for joining propagators 

In this section we consider a third way of calculating the equal space retarded propagator in 
the AdS-Vaidya background. Specifically we will consider the saddle point approximation, 
which is valid in the limit of large scaling dimension A, of the following integral formula 
for the CFT correlation function 



Gr(x 2 , t 2 ;xi, ti) — i 



dx dz 



99 



VfJ, 



v=0 



G^ s R (x } v 



0, z; x 1 ,t 1 ) 'd^ Gijf (x 2 , h; x, v 



0,z) 



(3.56) 



where t\ < < t 2 , and Gb,r and Gb,r are propagators that we will describe in more detail 
below. First, we have switched to coordinates (x, v, z) in which the AdS-Vaidya metric 
reads 



ds z 



1 



[- (1 - 6(v)R 2 z 2 ) dv 2 - Idvdz + dx 2 ] 



(3.57) 



The time coordinate t appearing in (3.56) is the boundary value of the Eddington-Finkelstein 
time v. 

In Lorentzian signature one has to distinguish between bulk-to-boundary propagators 
and boundary-to-bulk propagators. The retarded boundary-to-bulk propagator Gb,r hav- 
ing its initial point at the boundary and the final point in the bulk, evolves initial data from 
the boundary into the bulk. On the other hand, the retarded bulk-to-boundary propagator 
Gb,r takes bulk data and evolves it towards the spacetime boundary. Below we give a 
precise definition of these two-point functions in terms of the bulk-to-bulk propagator. 

In (3.56), the two sided derivative is understood to act only inside the square brackets, 
is the retarded boundary-to-bulk propagator in AdS 3 and is the retarded bulk- 

to-boundary propagator in BTZ. The integral in (3.56) is performed at the position of the 
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shell v — and g is either the BTZ or the AdS3 metric. 5 The basic logic is as follows. 
A delta-function source inserted on the boundary before the shell injection produces a 
disturbance that propagates to the shell according to Gg d R . After passing through the shell, 
the disturbance propagates back to the boundary according to G^ R Z . In the following we 
provide a derivation of (3.56). 



Derivation of Eq. 3.56: First consider the retarded bulk-to-bulk propagator Gbb,r 
associated to a scalar field of mass m in an asymptotically AdS spacetime. It solves the 
equations 

(□ - m 2 )G BB , R (x', v', z'\ x, v, z) = (□' - m 2 )G BBjR (x', v', z'\ x, v, z) 

= . 5{y -v')5(x-x')5(z- z') 

vanishes when v ' < v and satisfies the standard normalizable boundary condition at the 
AdS boundary. 6 

Given the initial data $>(x',v' = 0,z') = Q (x',z'), for v ' > we will show that the 
unique solution to the massive wave equation satisfying the boundary condition of vanishing 
non-normalizable mode is given by 

v', z') = i I dxdzyf ZI gg V11 ^q{x,z)oIg BBjR {x' ' -> v ' ,z';x,v = 0,z) . (3.59) 
Jv=o L 

For v' > 0, it is clear from the definition of the bulk-to-bulk propagator (3.58) that (3.59) 
solves the massive wave equation. By integrating (3.58) with respect to v from v' — e to 
v ' + e we get 

dv {d^{V--99^ v duG BB ^ R (x\ v', z'; x, v = 0, z)) - m 2 ^gG BBA {x\ v', z'\ x,v = 0, z)j 

= —i8(x — x')S(z — z'). (3.60) 

The equations of motion (3.58) imply that G BB) r is a regular function of v near v = v'. 
Thus, the mass term as well as the d x (y/^gg xx d x G BB ^Fi) term can be assumed to be regular 
near v = v ' leading to order e contributions to the integral in (3.60). The terms singular 
near v — v' in the integrand come from the derivative terms 



d v (V=gg vz d z G BB>R ) + d z (^g zv d v G BB , R ) } (3.61) 
which can be written as 

d v (V=99 vz d z G BB>R + d z {^gg zv G BBjR )) - d z {d v {^gg zv )G BBjR ) , (3.62) 



5 For the above formula to hold it is necessary for \J-gg" 11 to be continuous at the junction v = 0. As 
this is the case in the Vaidya spacetime, using either the BTZ or the AdS3 metric gives the same result. 

6 Here we focus on scalar fields satisfying the standard normalizable boundary conditions $ oc z A as 
z — >• 0, but generalizations to more general boundary conditions are in principle straightforward. 
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where we have dropped the explicit coordinate dependence. Using the fact that y/—gg vz 
is a continuous function of v, and assuming that d z G R is a regular function of v , 7 the last 
term above gives a contribution of order e to the integral in (3.60). Finally the integral 
(3.60) can be performed to give 



-gg vz d z G BB>R + d z (y/--gg zv G bb,r) 



-id(x — x')5(z — z'), 



(3.63) 



J v=v' —e 



where we used the fact that G R vanishes when v > v' . Thus, we see that the bulk-to-bulk 
propagator satisfies 

\ fzr gg v * 1 d^G BB)R (x ', v', z'\ x, v = 0, z) + d lJL (y/ zr gg v ^G BB)R (x' ', v', z'\ x, v = 0, z)) 

-> -i6{v')8{x-x')5{z- z'), (3.64) 

as v' — > 0; so that indeed (3.59) satisfies the correct initial condition <&(x,v = 0, z) = 
In particular we can use (3.59) to join two retarded bulk-to-bulk propagators 



G BB , R (x2,V2,z 2 ;x 1 ,v 1 ,z 1 ] 



i I dxdz^/-gg Vix 

>v=0 
^ y 

x d tl G BB}R (x 2 ,V2,z 2 ;x,v = 0,z 



G BByR {x, v = 0, z; xivx, , zi) x 



(3.65) 



To obtain the boundary theory retarded correlation function Gr from a bulk-to-bulk prop- 
agator we can use the well known relation (see e.g. [10,32]) 



G R (x 2 ,t 2 ; x 1 ,h) = lim (2u) {ziz 2 ) G BB , R {x 2} v 2 ,z 2 ;x u v u z 1 ) 
where t is the boundary value of v. From (3.65) we then obtain 



(3.66) 



G R (x 2 ,t 2 ;x 1 ,t 1 ) = % I dxdzy/^gg G B>R (x } v = 0, z;x 1 ,t 1 )d^G BiR (x 2 ,t 2 ;x,v = 0, z) , 

Jv=0 

(3.67) 

where we have identified the retarded boundary-to-bulk propagator as [ ] 



G B R (x, v — 0, z] X\,ti) = 2v lim z l A G BB R (x, v — 0, z\ X\,V\, Z\, 
and the retarded bulk-to-boundary propagator as 



(3.68) 



G B , R (x 2 ,t 2 ;x,v = 0,z) = 2u lim z 2 G BB! n(x 2 , v 2 , z 2 ; x, v = 0, z) . (3.69) 



7 This assumption is indeed consistent with the final result we obtain. 
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Application to AdS-Vaidya: In the specific case of AdS3 Vaidya, (3.67) reduces to 
(3.56). As reviewed in Appendix E, the retarded boundary-to-bulk propagator in AdS 3 
spacetime is given by 



z A 



G^(x, v = 0,z; x u h) oc 9(z - Vl )lm { — — s } , (3.70) 

[—t{ + 2zti + (x — Xi) 1 + ie\ 

and the retarded bulk-to-boundary propagator in BTZ is given by 



z A 



in terms of the coordinates 



G* T £(x 2 ,t 2 ;x,v = 0,z) cx 9(t 2 - z)lva\ A \ , (3.71) 

[ \-t z 2 + 2zt 2 + [x - x 2 ) 2 + ie\ J 

:dinates 

z = ze Rx , x = VR- 2 - z 2 e Rx cosh{Rt), i = VR~ 2 - z 2 e Rx sinh(itt), v = t-z. 

(3.72) 

The above coordinate transformation maps BTZ locally into AdS3, as reviewed in Appendix 
E. In the new coordinates, the shell is located on a surface v = where the coordinate 
relations reduce to x = ^ log Rx and z — zj (Rx). 

Substituting (3.70) and (3.71) and performing the change of coordinates (3.72), the 
integral (3.56) becomes 

Gr(x 2 , t2',Xi,ti) — J dxdze~ s , (3.73) 



where 
S = A 



2 log z + log [ Rxt\ — 2t\Z — Rx ( — log x — x\ 

+ log I -j^e 2Rx2 sinh 2 Rt 2 - 2^e X2 sinh Rt 2 - ( x - \e Rx2 cosh Rt 2 
R R \ R 



(3.74) 



to leading order in A. For simplicity here we consider non-integer A - the extension to 
integer A is described in Appendix C. The integration contour is determined by the ie 
factors in (3.70) and (3.71) and in general it can be rotated to pass through the saddle 
point 

d- z S = 0, d s S = 0. (3.75) 
For the equal space correlation function {x\ = x 2 = 0) the saddle point is at 

x = l/R, (3.76) 
2t 1 sinh(^) 



' Rtx cosh(^) + 2 sinh(^) ' 
Substituting the saddle point into (3.73) gives the retarded correlator 



(3.77) 



G fl (0;* 2 ;0,*i)cxe- s cx- : — \ : iDt _~ A , (3.78) 



[ti cosh(^) - |sinh( 
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Rt 2 

2 



which agrees (up to an undetermined normalization) with the geodesic approximation result 
(3.37) for t\ < < t 2 . To make this argument precise one would need to check that the 
saddle point is on a steepest descent contour, which is nontrivial in the presence of branch 
cuts. 

4 Beyond the geodesic approximation 

In this section we will develop an alternative method for computing propagators in AdS- 
Vaidya that goes beyond the geodesic approximation and enables us to compute the result 
for all spatial momenta k and for any conformal dimension A. This method allows us to 
compute the full spectral function but only gives numerical results. 

4.1 Holographic dictionary without time translation invariance 

Usually the holographic dictionary is discussed in translation invariant backgrounds. The 
Green functions then depend only on the difference x — x' of two points. When translation 
invariance is broken in at least one direction, the Green functions depend explicitly on two 
points. A standard holographic method to compute the retarded Green function in the 
field theory at the boundary is to study the dual bulk field: choose an infalling boundary 
condition at a (coordinate) horizon, solve the field profile from the equation of motion, and 
then extract the Green function from the asymptotic behavior. So we will begin with a 
brief discussion of the holographic dictionary for one-point and two-point functions in the 
case of broken translation invariance. 

Consider a free boson in 3d asymptotically AdS (AAdS) space of action 

S <P = -\J d 3 xV=9 {d^dy + m 2 2 ) + Set , (4.1) 

where the counterterm action is given by 

S c t = J dxdt^0 2 \ z=e , (4.2) 

and 7 is the induced metric on the z = e surface. 8 In the field theory on the boundary, 
using standard results of linear response theory applied to AdS / CFT [34-36] and large-N 
factorization, the response of an operator at (x 2 ,t 2 ) in the presence of a delta function 
source at (xi,ti) is given by the retarded correlation function 9 

(0(x 2 , t 2 )) s = i9(t 2 - ti)([0(a: 2> t 2 ), 0{x x , h)}) = -G R (x 2 , t 2 - xuh) . (4.3) 

8 The counterterm (4.2) is applicable for the range of scaling dimensions A < 2, to which we specialize in 
the following numerical calculations. For larger A, and for more complicated bulk field content, one needs 
to add more counterterms (see e.g. [33]). 

9 Here we are assuming that O does not have a vacuum expectation value. 
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The details of the derivation of (4.3) are included in Appendix F. On the other hand, 
according to the standard holographic dictionary, the one-point function in the boundary 
theory is given by 

£ ^fon— shell 

(0(x, t))j = + = -2u(f>+(x, t) , (4.4) 
oJ(x, t) 

where J and <p + denote respectively the leading and subleading asymptotic behavior of the 
bulk scalar field <fr as z — > 

<p(x, t, z) = J(x, t)z A - + <p+(x, t)z A + ... (4.5) 



2 



where A = 1 + u, A_ = A — 2v and v = y/l + m 

Comparing (4.3) and (4.4), we obtain the two point function 

G R {x 2 , t 2 ; x lt h) = 2v<j) +: s(x 2 , t 2 ) , (4.6) 

where <p + g is the subleading mode of a bulk solution that satisfies the boundary condition 
J(x2,t 2 ) = S(x 2 — xi)5(t 2 — t\) for the non-normalizable mode. A similar formula for 
the retarded correlator in momentum space in a time translationally invariant state was 
obtained in [36]. For more discussion on the relationship between the one point function 
procedure and the more conventional prescription for obtaining retarded correlators [35], 
see [36]. To construct such a solution, it is more convenient to specify the boundary data 
and reconstruct the full bulk field profile using the retarded boundary-to-bulk propagator 

(p(x 2 ,t 2 ,z) = J dxdtG^ A g s (x 2 ,t 2 ,z;x,t)J(x,t) = G^ s (x 2 ,t 2 ,z;x 1 ,t 1 ) , (4.7) 

rather than specifying the boundary condition at the Poincare or event horizon and inte- 
grating towards the boundary. 

In the case of AdS3 spacetime, the boundary-to-bulk propagator is given by 

G^t^M) = C^ h - t ^ Hx2 _ Xi) :_^_ ti)2 + i f , (4.8) 

where Ca = r(A)/(7iT(A — 1). Using (4.7) we obtain the solution 0, whose normalizable 
mode can be read off, which after using (4.6), leads to the CFT retarded two point function 

Gr uum (^t 2 ;x 1 ,t 1 )=g(t 2 -t 1 ) 2(A ~ 1)r( f ) Im( ? —— ) A , (4.9) 

7iT(A-l) \{x 2 - xi) 2 - [t 2 - ti) 2 + it/ 

This also gives us the two point function in the CFT state dual to the Vaidya spacetime 
before the shell for t\ < and t 2 < 0. 

In the next section we apply this method to the Vaidya background, the only difference 
with respect to the pure AdS setup being that to construct the bulk field profile one needs 
to continue the solution across the null shell. 
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4.2 Numerical construction of the CFT retarded two-point function 

As discussed above, to obtain the retarded boundary correlator, we need to solve the scalar 
equation of motion (□ — m 2 )cf) = in the Vaidya metric background. We parametrize the 
geometry using ingoing Eddington-Finkelstein (EF) coordinates and the metric (3.57) and 
define 

h(z,v) = 1 -6{v)R 2 z 2 . (4.10) 

Since the metric is translation invariant in the spatial coordinate direction x, it is 
convenient to perform a Fourier transform. The equation of motion for the scalar field 
then becomes 

hd 2 z <P + \hd v <f> - 2d v d z A + (d g h - -^j d z cj> - ( A ~ ( ^ 2 ~ 2) + k 2 ^j = 0. (4.11) 

This equation is first order in the EF time v, and the horizon is a characteristic surface. 
This means that a solution to the equations of motion is uniquely specified by a single initial 
condition at a constant v slice, together with a boundary condition at the AdS boundary. 
Also, this means that the only condition one has to impose at the position of the infalling 
shell v — is the continuity of the field <p. This makes the collapsing null shell simpler 
to analyze than the case of timelike collapse, where one needs to specify two independent 
junction conditions at the (time-dependent) position of the shell. 

Due to its causal nature, for t\ > the retarded Green function is automatically thermal. 
Here we focus on correlation functions in the boundary theory with the initial point (x±, ti) 
in the zero temperature (AdS vacuum) region (ti < 0). Note that it is simple to start with 
the initial condition in the vacuum region: the bulk scalar field solution before the collapse 
is given by the standard retarded AdS vacuum boundary-to-bulk propagator. It solves the 
equation of motion (4.11), it approaches a delta function at the boundary and vanishes 
outside the future lightcone of (xi,ti). 

The Fourier transform to mixed (k, v, z) space of the retarded AdS3 boundary-to-bulk 
propagator is (see Appendix G) 

G%%{h,v,z',t 1 )=4>{k,v,z) 

= -C 0("-*i)* A \ k \*+hj (\ kly / {v _ tl) 2 + 2(u _ h)z \ } (412 ) 

[{v-t l ) 2 + 2{v-t 1 )z}— 2V J 

where 



In particular, at the shell v = the field profile is given by 

<ftk,v = 0,z) = -C d thh_ r \ kr +hj (\k\Jtl-2hz) , (4.14) 

and (4.12) indeed satisfies the asymptotic behavior (4.5) with J(t) = S(t — ti): 

(f>{k,x,v) =5(t-t 1 )z A - + ■■■ , (4.15) 
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where we recall v reduces to t on the AdS boundary. 

In order to construct the full field profile everywhere, we must extend it into the region 
after the null shell collapse (v > 0). As we know the behavior of the field at v < 
analytically, the calculation reduces to solving the first order equation of motion (4.11) in 
the BTZ background v > 0, with initial data specified by (4.14). In the following we solve 
(4.11) numerically. For convenience we choose the radius of the event horizon to be R = 1 
by a conformal rescaling. 

For v > 0, the scalar field equation (4.11) can be solved noting that it can be written as 



d z 



d v (, 



2<Jz 



hdl6 + d z h 



d z 



A_(A_ -2) 



k 



which can be integrated with respect to z to give 



+ 



dz 



i\Tz 



hdU +[d 3 h 



did) - 



A_(A_ -2) 



k 



(4.16) 



(4.17) 



in terms of the regulated AdS boundary z = Zq. For the numerical results presented in 
Section 4.3, we have set z = 10~ 10 . 

For computational convenience we define a rescaled field P(k, v, z) as 



0(/c, v, z) = z A ~ P(k, z, v) 
in terms of which (4.17) becomes 



(4.18) 



d v P\ 



1-2A_ 
+ Z 2~ 



d v P\ 



dz 



Z 2 



2() 



hdjP +(d g h + 



(2A - l)h 



d z p - (a 2 + k 2 ) p 



To solve this equation numerically, we discretize both space and time 

(v,z) ->• {v[i},z[f\) = (ih v ,jh z ), 



(4.19) 



(4.20) 



where (h v ,h z ) = (T/N v , 1/N Z ), N V ,N Z are the numbers of lattice sites and T is the range 
of the time coordinate. The field P(v,z) is replaced by discrete variables 



P(v,z)^P[ l ,j} , 
derivatives are replaced by discrete derivatives 

R[i + l,j]-R[i-l,j] 



d z P(v, z) ->- 



2K 



and the integral by a sum 



J dz — >■ h z . 



(4.21) 



(4.22) 



(4.23) 
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By knowing the value of the field P at some initial constant v slice, we can use the dis- 
cretized version of (4.19) to obtain the value of the field at the next time step. Iterating 
this procedure we obtain the time evolution of the scalar field. Alternatively, a similar 
discretization and integration procedure is efficiently implemented by the Mathematica 
command NDSolve. The results of the two computations agree and the plots we present in 
the next Section have been produced using Mathematica. 

4.3 Numerical results 

4.3.1 Retarded two-point function at fixed momentum 

The left panel of Figure 7 shows the time evolution of the retarded two-point function 
for fixed momenta (in blue continuous), together with the time evolution of the vacuum 
(dashed red) and thermal (dot-dashed black) correlators. The two-point function in the 
thermalizing CFT decays faster in time than the vacuum correlator. 



\G R %t2,t x )\ 

l k"\- 




Figure 7: The left panel shows the correlator in the thermalizing background (solid blue 
curve), together with the vacuum (dashed red) and thermal correlator (dot-dashed black). 
The right panel compares logarithmic plots of the thermalizing correlator (solid blue) and 
of the exponential decay of perturbations close to equilibrium set by the imaginary part of 
the lowest quasinormal mode (orange dot-dashed). In both figures we have fixed A = 5/3, 
t\ = —2 and k = ir/2. 



A logarithmic plot reveals that the absolute value of the correlator oscillates, with 
the enveloping curve decaying exponentially in time (right panel of Fig. 7). Ingoing field 
fluctuations in a black hole background are well known to have a quasinormal behavior with 
a set of discrete oscillation frequencies with negative imaginary parts. These modes have 
been identified in the AdS/CFT correspondence as determining the approach to thermal 
equilibrium in the level of linear response [37]. Even though in our case we are considering 
dynamics far from equilibrium, it seems reasonable that at least at large times the dynamics 
of the field is dictated by the lowest quasinormal modes. The right hand side of Figure 7 
compares the exponential decay of the field modes with those of the lowest quasinormal 
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mode, which in the black brane background is known to be given by [35] 

Uquasi = k - iA. (4.24) 

Also the oscillation frequency of the correlator is given by the real part of the lowest 
quasinormal mode, here k. 

4.3.2 The momentum dependence of the retarded autocorrelator 

Numerical results for the thermalizing retarded correlation function as a function of mo- 
mentum, at fixed time, are shown in Figure 8. The plot reveals a modulation by the power 
law k u at large momentum which is also present in the conformal field theory vacuum 

GT uum (k,t 2] tt) oc -9(t 2 - *i) (t^t-) +2 J -u-^\k\it2 - ti)) , (4.25) 
which can be obtained by Fourier transforming (B.27). 
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Figure 8: Numerical result for the correlation function at a fixed time t% = 5 and t% = —2 
in the thermalizing (solid blue) and thermal (dot-dashed black) background, as a function of 
momentum. The vacuum curve is not shown in the plot because of much larger amplitude. 
In the figure we have chosen A = 5/3. 



For small momenta the thermalizing retarded correlator in Figure 8 differs significantly 
from the thermal correlator. In particular, the thermal correlator is a non- vanishing function 
of k at small k, while the thermalizing correlator vanishes at small k. This non thermal 
form at small momenta is in qualitative agreement with the "top-down" thermalization 
pattern obtained from many of the holographic thermalization studies. 

4.3.3 Spectral function 

As in Sections 2.2 and 3.1.3, we introduce a mean time T = t\ + 1 2 and a time separation 
t = t 2 — t\. Fourier transforming Gjt(k,T,u) with respect to t 

POO 

G R (k,T,u)= / dte^GRfctxh), (4.26) 
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we obtain the time-dependent spectral function 

p(k,T,u) = -2ImG R (k,T,u). (4.27) 

In Figure 9, we display the frequency dependence of the spectral function for fixed momen- 
tum and various average times T, obtained by Wigner transforming our numerical retarded 
correlator. Unlike in the harmonic oscillator example, the spectral function here is a smooth 
function with no quasiparticle peaks. In the vacuum the spectral function has the simple 
form (B.37): p oc 8(u — k)(u 2 — k 2 ) v . The momentum k determines the frequency when the 
spectral function obtains support. In the thermal state, the spectral function has support 
all the way down to u> = and at large u> approaches the vacuum form. The time dependent 
spectral function in Fig. 9 smoothly interpolates between the vacuum and thermal forms. 



p(k,T» p(k,T»-p thermal (k,w) 




Figure 9: Wigner transformed correlation function without (left) or with (right) the 
thermal two-point function subtracted. The different curves correspond to different values 
of the mean time T: the vacuum T = — oo (red, dashed), T = — 1 (orange, dot-dashed), 
T = (blue, solid), T = 1 (brown, dotted), the thermal T = oo (black, dot-dashed). In 
the figure we have set k = 1 and A = 11/8. Note that in empty AdS the spectral function 
vanishes for u < k as it should. 



5 Conclusions and outlook 

We have studied some aspects of thermalization following a holographic quench in the 
strongly coupled quantum field theory dual to the AdS 3 geometry In this scenario, energy 
is suddenly injected in the ultraviolet in the form of a thin shell of null energy that falls in the 
bulk AdS3 space eventually forming a BTZ black hole. In particular, we have investigated 
unequal-time two-point functions of the boundary field theory via their holographic duals, 
either in the semiclassical approximation by calculating the length of the bulk geodesic 
connecting the endpoints of the boundary correlator, or by solving the appropriate wave 
equation in the time-dependent bulk geometry. 
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One quantity of particular interest is the spectral function of the correlator, which 
changes from the vacuum spectral function to the thermal spectral function during ther- 
malization. While causality dictates that the vacuum spectral function p(u, k) exhibits a 
gap uvi n = k caused by the absence of space-like excitations, the thermal spectral function 
exhibits no gap due to the presence of real excitations in the thermal state. In the case of the 
quench followed by thermalization, the time-dependent spectral function smoothly changes 
from the vacuum spectral function to the thermal one. Because it involves a Fourier trans- 
form covering an infinite time range, the time-dependent spectral function deviates from 
the vacuum spectral function even before the quench, especially for low frequencies, where 
it does not exhibit a gap. 

The techniques developed in this article enable the calculation of many quantities of 
interest in the thermalization process. For example, one could study the time dependence 
of the occupation number of field modes in the boundary gauge theory. 
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A A notion of time- dependent entropy 

Here we summarize the arguments in [38] that lead to the entropy formula given in Sec. 2. 
Let us consider particles with Bose statistics indexed by a macroscopic quantum number 
/. Also suppose that the states index by / have a degeneracy gi, and that there are nj 
particles with quantum number /. In both quantum physics and in statistical physics we 
would really have a probability distribution over nj, but we will assume that nj ^> 1 and 
the fluctuations are much smaller than the mean so that we can ignore them. Thus we will 
take rti as given and will not integrate over the distribution of this variable. 

Then, the number of ways of distributing the nj particles over the gj microstates is 
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The numerator counts the number of ways of putting gj — 1 partitions in a set of nj 
particles and the denominator comes from the indistinguishability of the particles. Here we 
are assuming that the particles are bosons. Then the number of total possible states of the 
system is 

T({n I }) = l[W(n I ). (A.2) 

Now let a be a particular micro state with macro occupation numbers {nj}. Assuming 
that all the microstates are equally probable, 



r(K» 



(A.3) 



is the probability of the microstate a 

Thus we can write a formula for the von Neumann entropy 



(A.4) 



S = -^P W} (a) \nP {ni} (a)=\nT({nj}) = ^lnW(nz) 

: 1. Assuming that n I) g I > 1 we can use 

s = ^2( n i + 9i) ln ( n / + 9i) - 9i In 9i ~ n i In n u ( A -5) 



where we used the fact that ^ a Pj nj }(ci 
Stirling's formula to write 



where we also dropped the — 1 in ni + gi — 1 and gj — 1. We can simplify further by factoring 
out gi : 



There are some cancellations, and then we can write 



9i 



5 + + l 



9i 



9i J \ \9i 



m\ (ni 

in 



(A.6) 
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9i 



(A.7) 



Here ni/gi is the average occupation number of microstates. 

To make contact with the notions introduced in section 2, we interpret gi = p(k,u,T), 
i.e., we take gi to be the time dependent density of states with macroscopic quantum 
numbers k, to. We also take rij = n(k,u,T)p(k,u,T), i.e. the total occupation number of 
modes with k, u. We then get 



S(T) = / duo dkp(k, tu, T) s(k, u, T) 



(A.8) 



where 



s(k, u, T) = (n(k, u, T) + 1) In (n(k, u, T) + 1) - (n(k, w, T)) In (n(k, w, T)) . (A.9) 
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To illustrate this formula, consider the simple case when the spectral density is a sum 
of (on-shell) delta functions: 

p{k,u,T) = Y i 6{u-u> i {k)), (A.10) 

i 

where the sum may account for any degeneracy by having several Wj(k) have the same 
value. Then in equilibrium: 

n(k,u) = n B (u)p(k,u) = y]n B (cOiO<-))5(co - w;(k)), (All) 

i 

where = [eP^ — is the Bose distribution. Thus nj = ns{oj) in equilibrium only 

depends on u, not on k, and it is independent of the level density and the degeneracy of 
levels. The formula for the entropy then becomes: 

s(co) = (tib(w) + 1) ln(ns(w) + 1) — wb(w) lnn^u;). (A. 12) 

The quantity s accounts for the entropy associated with the different possible occupation 
numbers of one bosonic mode of energy u>. 



B Euclidean AdS 3 and BTZ 

B.l Geodesic length in Euclidean AdS 3 and BTZ 

We here explicitly derive the renormalised length of a geodesic between two boundary 
points in Euclidean BTZ and AdS3. This serves both to introduce the Wick rotation in 
this setting, as well as to derive a number of expressions that will be used in the main text. 
We parametrize Euclidean BTZ using Poincare coordinates in which the metric reads 

ds 2 = (r a _ R 2 )dy 2 + + r 2 dx 2^ (R1) 

where R < r, the surface r = R is the black hole horizon and r — oo represents the 
spacetime boundary. Since Euclidean AdS3 is recovered from (B.l) by sending R — > 0, we 
will work out explicitly the BTZ case alone and take the AdS limit only at the end. 
The Euclidean BTZ geodesies are given by 

r (A) 2 = T 2 +(r^-r 2 )cosh 2 (A-A ), (B.2) 

1 /T_ \ 

x(X) = Xq + — arctanh I — tanh(A — A ) , (B-3) 
R \T + J 



2/(A) = y + -^arctan I A / p 2 _ ^ 2 tanh(A - A ) [ , (B.4) 
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for constants < T 2 < R 2 < and Xq, yo, Xq. As a function of r, we find the following 
branches (r G [r + ,oo[), 



1 /r_ r 2 -Y 2 + , 

x±{r) = x ± — arctanh I — J ^ _ ^ 2 J , (B.5) 



1 / hp _ y 2 r 2 - T 2 

y± (r) = y ± -arctan J p2 _ R2 J r2 _ T t I > < B ■ G > 



A ± (r) = A ±ln| ^ J^ff^ 1- (B.7) 

The separation between two boundary insertion points (xi,yi,r = oo) and (ar 2 ,2/ 2 ,r = oo) 
is given by 

2 /r \ 

Ax = x(X oo) — x(X — oo) = —arctanh I — , (B.8) 

r \ r i / 



2 / /i? 2 — r 2 

Ay = y(X ^ oo) - y(A -> -oo) = -arctan I J r2 _ ^~ ) , (B.9) 



and the length of the geodesic can be found to be 



2 



AA A + (r-M-„)-A-(/--M-„) = 2h.[ v u ^ 2 _ r2 " ± j (B.10) 

in terms of a regularized AdS boundary at r = ro- In the limit ro — >• oo this reduces to 

AA = 2 In ( ; 2r ° 1 . (B.ll) 

We now pause for a moment and consider the R — > AdS limit. In this limit, the 
geodesies branches (B.5)-(B.7) become 



1 / E 2 

x± (r) = x ± -cos(cr)Wl - — , (B.12) 



y±(r) = y ±-sm(a) ] Jl-—, (B.13) 

A±(r) = A ± arccosh (£) = A ± In ^v^-E 2 ^ (R M) 

where r G [E, oo[ and 

r 2 

lim T 2 , = E 2 and lim — - = cos 2 (cr) . (B.15) 
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The geodesic length (B.ll) thus reduces for R — > to 

AA = 2 In (^j = 2 In (WtAx) 2 + (Ay)*) . (B.16) 

We renormalize the geodesic length by subtracting the divergent contribution 21nr , so 
that the renormalised length of the geodesic in Euclidean AdS is given by 

6C EAdS = In [(Ax) 2 + (Ay) 2 } . (B.17) 

In a similar way, we renormalize the Euclidean BTZ result (B.ll) and obtain 

5 (^)w(^); 



<5£ EB tz = AA - 21n(r ) = In ( ^ 2 _ p2 ) = In 



(B.18) 



B.2 Two point functions of vacuum and thermal CFT 

The Euclidean two-point functions can be computed in the geodesic approximation from 
the renormalised geodesic lengths that were obtained in the previous section. A suitable 
Wick rotation leads to the time-ordered and retarded two-point functions in Lorentzian 
signature. Below we also briefly comment on how to extract the spectral function and 
occupation number from these results. 



B.2.1 Euclidean two point functions 

Remember that, in the geodesic approximation, the renormalised two-point function in the 
dual CFT equals 

(0(x u yi )0(x 2 , y 2 )) rcn ~ e- ASC . (B.19) 
From (B.17) and (B.18), we then find the (Euclidean) vacuum two-point function 

^vacuum/ X = _ (B.20) 

E 1 ' {(Ax) 2 + (Ay)*] A 1 ' 

and the (Euclidean) thermal two point function 

^thermal, x = 1 /g ^ 

^(smh 2 (^)+sin 2 (^))] A 

where the first can be obtained from the second in the R — > limit of vanishing black 
hole mass. Despite the fact that we have used a geodesic approximation, the result for 
the vacuum and thermal two-point function is exact, being fully constrained by conformal 
invariance (apart from an overall scaling). 
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B.2.2 Lorentzian two point functions 

To go from Euclidean coordinates (x,y,r) to Lorentzian coordinates (x,t,r), we perform 
the Wick rotation y = it. It is well known that by Wick rotating the Euclidean two point 
function, one finds the time-ordered (Feynman) two point function 

iG F (x 2l t 2 ; xi, ti) = G E (x 2 , it 2 ; x 1: it x ) . (B.22) 
For non integer scaling dimension A, this leads to 10 
iG—(x fx t) - ^) 2 -(^) 2 l , dmy-(Axf] 

" [(Ax)2 _ (At)2] A + [(At)2 _ (Ax)f e > 

*r"fx fx t) - ^[(Ax) 2 - (At) 2 ] 



[A (sinh 2 (^)-sinh 2 (^))] A 
[(At) 2 - (Ax) 2 ] 



e" (B 24) 

' ^(sinh 2 (^)-sinh 2 (^))] A 

The identity, 

G F (x 2 ,t 2 ] x 1} tx) + (G F (x 2 , f, xx,tx))* = G R (x 2 , t 2 ; xx, h) + G A (x 2 , t 2 ; xx,tx) , (B.25) 
where Gr/a denotes the retarded (advanced) two-point function, allows to deduce 

G R (x 2 ,t 2 ;x 1 ,t 1 ) = 9(tx - t 2 )[G F (x 2 ,t 2 ;x 1 ,t 1 ) + {G F (x 2 , t 2 ; x u tx))*) . (B.26) 
In our case, this results in 



9 [(At) 2 - (Ax) 21 
[(At) 2 - (Ax) 2 ] 



GT num (x2,h;xx,tx) = -2sin(nA)9(At) r/ [ lJ - \ J > (B.27) 



^thermal/ f f ) = _2 sin (ft A) 9 (At) jr. (B.28) 

J ^(sinh 2 (^)-sinh 2 (^))] A 

B.2.3 Spectral function 

Performing a Fourier transform 

/+oo /"+00 

dx / dte ikx e iult G(x,t), (B.29) 
-oo J —oo 

where we take x = Ax = x 2 — X\ and t — At — t 2 — tx, we can calculate the spectral 
function p(k,u) and occupation number n(k,u) from the identities (2.3) and (2.4) 

p(k,co) = -2ImG R {k,u), (B.30) 

(1 + 2n(k,u))p(k,uj) = —2 lmG F (k, u) . (B.31) 
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For the case of integer A, see Appendix C. 
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The Fourier transforms of (B.23), (B.27), (B.24) and (B.28) read respectively 



G vacuum 



(k,u) 



r(i - A) 
r(A) 



u 2 - k' 



2\ A-l 



9 (oo 2 - k 2 ) [isin(TrA) - cos(ttA)] 



, 2 o\ A-l 

k z — u zx 



6(k 2 -u 2 ) 



^B.32) 



G vacuum 
R 



(k,u) 



r(i - A) 
r(A) 



u 2 - k' 



2\ A-l 



9 (cu 2 — k 2 ) [i sign(o;) sin(7rA) — cos(7rA)] 



+ 



7 2 2\ A-l 

k z — u zx 



e(k 2 -uj 2 ) 



(B.33) 



G thermal 
F 



R 



2A-2 



2sin(7rA)(r(A)) 2 



r 



cos(7rA) cosh y~^r j ~ cosn 



G 



thermal 
R 



(k,u) 



R 



2A-2 



2sin(7rA)(r(A)) 2 



r 



A i(u + k) 
~2 + 2R 

nk 
~R 

A i i(u + k) 
2 



A i(co - k) 
~2 + 2R 



% sin(7rA) cosh 



/TTUJ 



2R 



r 



V R 

A t{u-ky 2 
~2 + 



(B.34) 



'7TUA 



cos(7rA) cosh ( 

V R J 



) — cosh 



nk 
~R 



2R 

— i sin(7rA) sinh 



(B.35) 



and agree with those computed in [35] up to an overall normalization. 
From these expressions, we can deduce that 11 

A-l 

Pvacuum 



Pthermal 



(T(A)) 2 
R 2A ~ 2 



4 



r 



A i(u + k) 
~2 + 2R 



(r(A)) 2 

Integrating over all momenta gives the results 

4vr 2 



9 [uj 2 — k 2 ) sign(o; 

r 



A i(u - k) 
~2 + 2R 



sinh 



Pv 



2A-1 



Pthermal( w ) 



T(2A) 

4?r J R 2A - 1 

r(2A) 



sign J 



^B.37) 



(B.38) 



r a + 



~R 



Sinh (if) 



(B.39) 



while for the occupation number, using (B.31), we recover Bose-Einstein distribution 



^thermal 



1 /ImG 



thermal 



(k,u) 



2 VImG^ ermal (fc,u;) 



2ttoj 

e r — 1 



L For half integer values of A, we can use the following expressions: 

2 



|r(l + ia:) 



sinh(7ra:) 



1 



IX 



cosh(7T2:) 



and T(x + 1) = xT(x). 



(B.40) 



(B.36) 
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C Notes on the case of integer Delta 
C.l Euclidean signature 

Consider a conformal field theory in d dimensions. The Euclidean two-point function of a 
scalar operator with conformal weight A is given by 

The Fourier transform of this function is 

G(k) = J d d xG(x)e !k ' x = J d d xe !kx | dss^e' 8 ^ 

1 f°° , a_i /7T\I _^ vri 



i A-l / " \ 2 -s- " / i a-A-l - — 

dss - e 4 S — — — / d«!i 2 e < 



r(A)i v s ; r(A)i 



r(A) 

where we have used the Gaussian integral 



/ 



dxe ikx e-P x2 = A /|e"fc , (C.3) 



the relation (valid for a > 0) 

— = W/ — 7 / dss^e'^, (C.4) 

and made the substitution s — 1/u 

This derivation is only strictly valid for d/2 > A > 0, which is not of physical interest 
to us since we always have A > d/2. However, one can analytically continue the result in 
A to all positive real values for which A — d/2 is not zero or a positive integer. We could 
also use a regulator £, such that 

<? € (x) = (x2+ \ r (C5) 



and determine the Fourier transform 

G € (k) = | d d x^(x)e ikx = -L^ | d^e^J^dss^ 1 



e -^ 2 e -< 2 



(A)/ 



OO do _^ 



ds s A 1 - e 4, e sl; = —— / — — -, — e * u e « 



27T5 / v 7 ^ 



T(A) V 2£ / 



r(A)7 u a-i+i 



K A _ d (VtfO, (C6) 
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where we have used the following expression for the modified Bessel function (valid for all 
a and b) 



K v (2ab) 



ds 

— I exp ( -sa 



(C.7) 



If v = A — d/2 is non- integer, then we can use the following series expansion around x — 0. 

, 2k 



KJx) 



7T (I- V (x) ~ Iy(x)) 

2 sinfz/Tr) 



and 



Thus we find that, 



IJx) 



2 



E 

k=0 



(f)' 



kW(u + k + l)' 



G c (k) 



2tt5 
f(A) \~2l 



^(vk 2 e) 



7T2 



T (A) sin (i/tt) 



E 

fc=0 



k 2 g 2 

4 



fc!r(-i/ + jfe + 1) 



2 \ oo 



E 

k=0 



k 2 g 2 

4 



A;!r(^ + A; + 1^ 



(C.9) 



If z/ < 0, then the limit £ — > of this result converges towards 



7T2 



+ 1 



G(k) = lim/ £ (k) = -r N , 

v; ^o m; r(A)T(z/+ l)sin(z/vr) V 4 



7T^ 



-r(f-A) /kV 



r(A) 



. (CIO) 



This correspond exactly to the result in (C.2). If v > 0, the first terms will diverge in the 
limit £ 0. However, note that these terms are analytic in k 2 ; they correspond to contact 
terms and are of no physical relevance. The only non-analytic term that survives in the 
limit £ —7- is given by 



G(k) 



7T2 



+ 1 



r (A) T(u + 1) sin(W) \ 4 



2\ v 



7T 2 



# r (f-A) /k^ A - 
r(A) \a 



(C.ll) 



This is the same result that we obtained by analytic continuation of (C.2) in A. If n 
A — d/2 is a positive integer, then we can use the series expansion 



K n {x) 



1 fx 

2 V2 

1 



n-1 

E 

fc=0 

'DO 



(n-Jfe-1)! 
k\ 



x 

T 



2\ fc 



+ (-l) n+1 In 



-^n(^) 



(ij>(k + 1) + i>(n + k + 1)) fx 2 ^ 1 



k=0 



k\{n + k)\ 



(C.12) 
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to write (using also the series expansion of I n (x) as above), 
G € (k) = 




4 / ^ fc!r(ra + fc + l) 

+(-t) E^s^i^f ) >■ (cud 



fc=0 

k 2\™ ~ (^(A; + 1)+V(n + A; + 1)) /k 2 f x *' 



A;!(n + A;)! 

In the limit £ — > 0, the first terms will diverge; however, these are analytic in k 2 , so we 
can neglect them. The last terms are also analytic in k 2 . The only non- analytic term that 
survives in the limit £ — > is given by 

c (k) = a^! (h!Y" lll(k2) = i^±g , (?fW) • (c.i4) 

yj r(A)n! V 4 / r(A)r(A-| + i) V 4 / 

There is yet another way to find this last result. If A is such that n = A — d/2 is an integer, 
we could also have used another regulator 5 and replace everywhere A by A — 5. Then n 
gets replaced by n — 5, which is no longer integer. Now we can use that 

T (-n + 5) = (1 + ^{n + 1)5 + O (5 2 )) , (C.15) 

n\ 5 

' 'l + ii(A)6 + 0(6 2 ),) (C.16) 



r (A - s) r (A) 

,2\ n-S /, 2 \ " / /i,2 



(t) = (t) (i-fc(£)' + <>(0), (C17) 



such that 

G(k) = 7T 



,r(f-A + 5) /k 2 \ A ~^ 



r (A - s) v 4 



' 1 ' • c(n • 1) • 0(A) In , ) • 0(5)). (C.18) 



r(A)n! V 4 / \$ 

The term in 1/5 is divergent and analytic. All terms of order 0(5) vanish in the limit 
5—t-O. The non-divergent term contains an analytic piece and a non-analytic one, namely 

*^w(T)" ,l(k '»' (C19) 

which is exactly the right result. 
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C.2 Minkowski signature 

We can perform a Wick rotation of the Euclidean two-point function 

GW = (x^ = (j,2 +£2) A' ( C - 20 ) 

to Minkoswki signature by taking y = e %e t and continuously take 9 from to n/2 — e (where 
e > and e — > is implied). We then find the time-ordered two-point function 

«W) = pFTiW (a21) 

For non-integer A, we can rewrite this as 

1 9 (x 2 - t 2 ) 9 (t 2 - x 2 ) 



iG F (x) = 



£2_£2)A (f2_ t 2)A (_( t 2_£2))A 

ITT A 



(x 2 - t 2 ) 9 (t 2 - x 2 



(C.22) 



(£2_ t 2)A ( t 2_ f 2)A- 

For integer A, one needs to pay attention to contact terms that may arise as for example 
in 

«Mx) = p_) 2 + ie = - in5 if - t 2 ) . (C.23) 

In the mean text of this paper, we have assumed that A is non-integer. One can find the 
results of integer A by using the regularisation A — 5 as in the previous section, where 
the contact terms result in analytic terms in k 2 . One can show that this will change 
the expressions for the real part of the time-ordered and retarded propagators, but not 
the imaginairy part. Therefore, the expressions for the spectral function and occupation 
number are valid for all values of A. 

D Complexified geodesies in AdS3 and BTZ 

In this Appendix, we review the construction of complexified geodesies in AdS3 and BTZ. 
AdS 3 

In empty AdS3 spacetime 

ds 2 = -r 2 dt 2 + ^+r 2 dx 2 , (D.l) 



the (real) geodesies are given by 
• Spacelike geodesies 
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* j 2 > e 2 : 



V(A) = VF^cosh(A-Ao) 

x(X) = x + ^ tanh(A - A ) (D.2) 
t(X) =t + tanh(A - A ) 



-k e 2 > j 2 : 




j 2 sinh(A - A ) 

^coth(A-Ao) (D.3) 
ji coth(A - A ) 



Timelike geodesies 



-k j 2 > e : no solutions 



-k e 2 > j 2 : 



r(A) = \/e 2 — j 2 cos(A) 

x(A)=x + ^tan(A) . (D.4) 



t(X) = t + ^zj2 tan(A) 



Timelike geodesies oscillate in the bulk and do not extend to the boundary. The first 
class of spacelike geodesies with j 2 > e 2 go from the boundary into the bulk, and back to 
the boundary connecting spacelike separated points on the boundary. The second class of 
spacelike geodesies with e 2 > j 2 go from the boundary (r — > oo) to the Poincare horizon 
(r = 0). These are the ones that can be used to connect timelike separated points on the 
boundary. However, the expressions x(X) and t(X) diverge when the geodesic reaches the 
Poincare horizon. To regularize this behavior, we add a small imaginary part to the afline 
parameter A, so that the coordinates can thus in general take complex values: 

'r(A) = v /e 2 -j 2 sinh(A - A + 

x(X) = x - ^i-z coth(A - A + i/3) . (D.5) 
t(X) = t - ^ coth(A - A + if3) 

These expressions can also be rewritten as 

'r(A) = y/e 2 - j 2 (sinh(A - A ) cos(/3) + i cosh(A - A ) sin(/3)) 

j ( sinh(2(A-A (1 ))~isin(2/3A 

X{/\) — X e 2_p \^ CO sh(2(A-Ao))-cos(2 / 3) / , { U -'°) 

+ (\\ _ 4- _ e ( sinh(2(A-Ao))-isin(2^) \ 
, V'V L e 2 -j 2 \ cosh(2(A-A ))-cos(2 / 3) J 

where we see that in order to have x(X) and t(X) finite everywhere, we need cos(2/3) ^ 1. 
Then: 

Ax = lim x(X) - lim x(X) = — (D.7) 

A^oo A^-oo e 2 — J 2 

2e 

At = lim t(X) - lim t(X) = — = (D.8) 

A^oo A^-oo e 2 — J 2 
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such that At 2 — Ax 2 = 

e 3 

Note that for complexified r-coordinate, going from the boundary (at infinity) to the 
boundary (at infinity) corresponds to |r(A)| 2 going from the point at infinity in the complex 
plane back to the point at infinity. The direction in which this happens is determined by 
the parameter (3. We now rewrite the affine parameter A as a function of |r|: 



A±(|r|) = Ao ± arccosh 
so that the geodesic length can be expressed as 

AA(|r|)=A + (|r|)-A_(|r|) = 21n 



e 2 — j 2 



+ cos 2 (/3) 



(D.9) 



e 2 — j 2 



+ cos 2 



+ 



e 2 — j 2 



-sin 2 (/3) . (D.10) 



After subtracting the divergent part, we find the renormalised geodesic length 

4 



5£ Ad s=lim (AA(|r|)-21n(|r|)) = ln , 

|r|^oo \e z — J 



In (At 2 - Ax 2 ) 



(D.11) 



This result is independent of the regularisation parameter (3, and it agrees with our expec- 
tation for the scaling of two-point functions in a vacuum CFT. 



BTZ 



We can repeat the construction in BTZ spacetime 



ds 2 = 



R 2 )dt 2 + 



dr 2 



+ r dx . 



r 2 -R 2 

For spacelike geodesies, we now distinguish the following cases: 
* (J 2 - E 2 + R 2 ) 2 > (2JR) 2 : 

>(A) = \ ((J 2 - E 2 + R 2 ) + cosh(2(A - A )) y/(J 2 - E 2 + R 2 ) 2 - 4J 2 i? 2 ) 

, x{\) = x + ^arccoth [j^g^ - yj (*=gg&) 2 - l) tanh(A - A )) 

t(\) = t + ^arccoth ^^g=^ - J2 -^ R R2 ) 2 ^ tan h(A - A 



(D.12) 



(D.13) 



* (2Ji?) 2 > (J 2 - E 2 + i? 2 ) 2 : 

V 2 (A) = | ((J 2 - £ 2 + R 2 ) + sinh(2(A - A )) ^4 J 2 R 2 - (J 2 - E 2 + R 2 ) 2 } 
< x{\) =x + larccoth [^g^ tanh(A - A ) - yjl - (^g^)^ 

f (A) = t + ^arccoth (^g^ tanh(A - A ) - y 1 ~ (^§T^) 2 ) 



(D.14) 
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Note that arccoth(rr) = §hi(^) is only denned for x G K\] — 1, 1[. Nevertheless, we 
can take the analytic continuation arccoth(2;) for z G C\{ — 1, 1} which has a branch cut 
between the points 1 and -1. 

To be able to take the AdS limit R — > 0, we have to consider the first class of solutions 
with (J 2 — E 2 + R 2 ) 2 > (2JR) 2 . In this first class, we can adopt the convenient notation 



A± = X - ((J 2 - E 2 + R 2 ) ± a/ (J 2 - E 2 + i? 2 ) 2 - 4J 2 /? 2 ) , 



(D.15) 



such that A + A_ = J 2 R 2 and (A+ - i? 2 )(A_ - R 2 ) = E 2 R 2 . We then find the complex 
solutions 

>2(A) = A_ + (A+ - A_) cosh 2 (A - A + i(3) 

x(X) = xq + -^arccoth ^^/^tanh(A — A + i(3)j . (D.16) 
t(A) = t + ^arccoth ^ ^~Z% tanh(A - A + 
The separation on the boundary is given by 



and At = — arccoth 
R 





- R 2 \ 




-R 2 j 



which gives 



We also have 



which results into 



A_ 
A + 

r 2 - A_ 



R*COSh 2 (K^) 



sinh 



2 / RAx 



R 2 sin h 2 (^) 



A, - A. 



A-t = A ± arccosh 



sinh 2 - S inh 2 • 

cosh 2 (A- A ) -sin 2 (/3), 



r 2 - A_ 



A^ - A_ 



+ sin 2 (/3) 



and 



5£btz = lim (AA(|r|)-21n(|r|))=ln 

|r|— >oo \ | y V_4_ — A_ 



4 



= In 



R 2 V Sinh2 (^T 



sinh^ 



RAx 



(D.17) 



(D.18) 
(D.19) 

(D.20) 
(D.21) 



(D.22) 



as expected. 
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E A useful coordinate transformation of BTZ 

The metric of a planar BTZ black hole is given by 



ds 2 



-(1-R 2 z 2 )dt 2 + 



dz< 



1 - R 2 z 2 



+ dx z 



(E.l) 



It is locally identical to the metric of AdS3 as can be seen by performing the change of 
variables 



z = ze 

The relations (E.2) can be inverted to give 



x = VR- 2 - z 2 e Rx cosh(ift) , 
i = VR- 2 - z 2 e Rx sinh(ftt) , 

Rx 



(E.2) 



X= 2R hg [ R2 (z 2 + x 2 -? 

1 fx + i 

t — — log = 

2R & \x-t 



z 



(E.3) 



RVz 2 + x 2 -t 2 
and the BTZ metric (E.l) becomes 



ds 2 



1 



-dP + dz 2 + dx 2 



(E.4) 



In particular, the retarded bulk-to-bulk propagator in the BTZ black hole is simply 
identical to the vacuum one in the new coordinates 



—A W A+1 A -A-^2 

o 2 F 1 [^—,—,A,a 



iGl T B z R {x, t, z; x', F, z') = A6(t - t')lm 

where A is a real constant depending on A and we have defined 

\(z' z (x - x 1 ) 2 - (t - t'f + ie 



a = -I — + 

2 \z z' 



zz 



(E.5) 



(E.6) 



Using the definitions (3.68)-(3.69) given in Section 3.3, we obtain the boundary-to-bulk 
propagator 



iGl T g( x,t,z;x', F) = 2 v+2 uA6{t - t')lm j 
and the bulk-to-boundary propagator 
iGl T £{x, i; x, f, z) = 2 u+2 uA6(t - i')lm 



[-(t - t') 2 + {x- x') 2 + z 2 + ie] 



A f ' 



(E.7) 



(z') A 



■it - t') 2 + (x - x') 2 + (z') 2 + ie} 1 



A ■ (E.8) 
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F Relating Gr to the one-point function 

Consider a quantum system described by the following Hamiltonian 

H = H -J(t)O, (F.l) 

where we have suppressed the spatial coordinates as they play no crucial role in the follow- 
ing analysis. For concreteness we work in the Schrodinger representation where the only 
operator depending on time is the density matrix. In the following we will take H to be a 
generic Hamiltonian which is allowed to have explicit time dependence. 

Let us assume we have solved for the density matrix p when J = 0, meaning that we 
know the general solution to the equation 

id tPo (t) = [H ,po(t)]. (F.2) 

The solution po(t) can be written in terms of a time evolution operator U, which is defined 
by 

id t U{t) = H U(t), (F.3) 

so that 

p (t) = U(t)p (t = 0)rf(t). (FA) 

The idea of linear response is to work out the time dependence of the density operator to 
linear order in J(t). Consider the ansatz p = p + Sp, where we assume Sp to be of the 
order of J. The equation of motion for Sp at linear order becomes 

id t 8p(t) = [H , 8p(t)} - J(t) [O, p (t)} . (F.5) 

The solution can be obtained through the method of variation of constant 5p = UsW, 
which allows to integrate the equation 

s(t) = s(t ) + i f dt'Jit') [0 H (t'),po(t')\ , (F.6) 

J t 

where we have defined the Heisenberg picture operator Oh if) = U(t)OW(t). Assuming 
J(t) vanishes for t <t for some to we get the result for the density matrix 

p(t) = p (t)+t f dt'J(t')U(t)[O H (t'),p (t')]U\t). (F.7) 

J — oo 

The expectation value of the operator O becomes 

(0(t))j = Tr(p(t)0) =Tr(p (t)O)+i J dt' J(t')Tr (p (t) [0 H (t), H {t')\ ) , (F.8) 
so that the overall effect of the source J on the expectation value of O is 

/oo 
dt'J(t')G R (t,t'). (F.9) 
■oo 
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For the specific source J(t) = eS(t — to), one obtains 

(0(t)) eS - (O(t)) = -eG R (t, t ). (F.10) 
Dividing by e and taking the limit e — > gives 

^(O(t)) eS = -G R {t,t ). (F.ll) 

On the other hand, in the holographic setup the one-point function in the presence of a 
source of the form J = e5(t — to)S(x — x ) is given by 

X Con— shell 

where = e5(t — t )S(x — :r ),2 A - + z A +(f) + + ... as z — > 0. The dependence of <p + on e is 
linear since the bulk equation of motion is linear (in the leading large- iV approximation). 
Thus, 0+/e is independent of e and we have 

^-(0(x,t)) tS = 2u^±. (F.13) 
oe e 



Since the final result is independent of e we can set e = 1, which gives 

= {0(x,t)) s , (F.14) 



d 

G R (x,t;x ,t ) = —{0(x,t)) eS 



e=l 



where we have reinstated the spatial dependence. 

Notice that the assumption that <j) + is a linear function of e can fail if the operator O 
has a non-vanishing expectation value. We will not consider such a case further in this 
paper. 

G AdS3 retarded bulk-to-boundary propagator 

The Euclidean AdS boundary-to-bulk propagator is given by 

z A 

G^ d l(x,y,z;x',y') = C A x , (G.l) 

where 

c ^^Bt } - < a2 > 

This can be Fourier transformed into a mixed momentum space as 

/oo A p —ikx 
-oo {z + {y -y ) +x 2 )* 
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One way to perform this integral is to use the Schwinger parametrization 

/T POC POO 

G%%k,y,z-y') = ^ / dx ds s *-i z * e ^-<* 2 Hv-v') 2 +* 2 ) ? ( G .4) 
and perform the Gaussian integral over x to get 

fi ~A poo 2 

G A *l(k, y, z- y') = / ds s^e'^ , (G.5) 



+ Jo 

where we have introduced a = z 2 + (y — y') 2 . Introducing a new integration variable u = as 

we get 

G^i(k,y,z;y') = ^^a^l duu^e'^ , (G.6) 
which can be recognized as an integral representation of a Bessel function 



G^(k,y,z;y') = - ^k^a-^^K^k^ ■ (G.7) 

Next we want to compute the retarded propagator, which can be obtained from the Feynman 
one through the identity 

G BtR (k, t, z; t') = -9(t - t') (G B , F (k, t, z; f) + G* B F (k, t, z; f )) . (G.8) 

To determine the Feynman propagator in the future lightcone, we should continue (G.7) to 
real time as a = z 2 + (y — y') 2 — > —(t — t') 2 + z 2 + ie. Since we need fractional powers of a we 
need to take into account the ie factor to determine the phase a a -> {-(t-t') 2 + z 2 + ie) a = 
e lna b, where b is the positive quantity b = (t — t') 2 — z 2 . In continuing to real time the 
Feynman propagator is multiplied by a factor of i as a part of its definition. In this way 
we obtain 

G^(k,t,z;t') = ^z 1+ "\krh-^^H^ .(IfclVft), (G.9) 

Z " 2 

where we have defined a new normalization constant 

The retarded bulk-to-boundary propagator can be obtained from the Feynman propagator 
through the identity 

G%$(k,t,z;l!) = -C6(t -t' - z)z 1+u \k\ u+ h- 2 -^ J_ u _ i(Vb\k\), (G.ll) 

where the 6{t — t' — z) appears because G B d p is purely imaginary for b < 0, which is outside 
the future lightcone. 
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